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ABSTRACT 

The heat transfer behavior of accreting ice 
surfaces in natural (flight test) and simulated 
(wind tunnel) cloud icing conditions have been 
studied. Observations of wet and dry ice growth 
regimes as measured by ultrasonic pulse-echo 
techniques were made. Observed wet and dry ice 
growth regimes at the stagnation point of a 
cylinder were compared with those predicted using 
a quasi steady-state heat balance model. A series 
of heat transfer coefficients were employed by the 
model to infer the local heat transfer behavior of 
the actual ice surfaces. The heat transfer in the 
stagnation region was generally inferred to be 
higher in wind tunnel icing tests than in natural, 
flight, icing conditions. 

NOMENCLATURE 

A,B experimentally derived constants (-) 

C. specific heat capacity of ice (J/kg K) 

C^ specific heat capacity of air (J/kg K) 

specific heat capacity of water (J/kg K) 

D w Diffusion,coeff icient of water vapor in 

air (m/^/s ) 

d cylinder diameter (m) 

h local convective heat transfer coefficient 

(W/m z K) 

k thermal conductivity of air (W/m K) 

latent heat of fusion of water (J/kg) 

L $ latent heat of sublimation of water (J/kg) 

L s latent heat of vaporization of water 

(J/kg) 

M" local mass flux/time (kg/m s) 

Nu Nusselt number (-) ? 

Q" local heat flux/time (W/nr) 

Re Reynolds number based on cylinder d i a . 

and V w (-) 

r recovery factor, 0.875 (-) 

AT^ cloud supercooling = -T a ("C) 

T f equilibrium surface temperature ("C) 
r u cloud temperature (*C) 

t" icing time (s) 

freestream velocity (m/s) , 

W°° cloud liquid water content (g/rrr) 

B local collection efficiency (-) 

c*, saturated, vapor density over surface 

'• 5urf (kg/m 3 ) 

o saturated, vapor density in cloud 

v ’ (kg/rrr) 
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INTRODUCTION 

Whenever an aircraft encounters liquid 
moisture in the form of supercooled cloud 
droplets, or freezing rain, ice will form on the 
exposed surfaces. Typical cloud droplet diameters 
range from as large as 50 microns to less than 10 
microns. In the case of freezing rain, droplets 
may be several mm in diameter. The shape of the 
accreted ice and its affect on the aircraft's 
aerodynamic performance 1 depend on several 
parameters - the cloud temperature, the average 
cloud droplet size and size spectrum, the amount 
of liquid water per unit volume, W, contained in 
the cloud and the size, shape and speed of the 
accreting surface. 

RIME ICE 

“Cold" cloud temps. M0*C to -30*C) 
Small droplets (12p typ.) 



GLAZE ICE 

“Warm" cloud temps. (0 # C to -10*C) 
Large droplets (20p typ.) 



Figure 1. Typical "rime" and "glaze" ice shapes. 

When all of the impinging droplets freeze on 
impact with the accreting surface the ice 
accretion is character ized as rime ice. The 
resulting ice shape typically protrudes forward 
into the airstream (see figure 1). Relatively 
cold cloud temperatures (below -10*C) and small 
droplet sizes promote rime ice formations. At 
warmer cloud temperatures and for the 
characteristical ly larger droplet sizes and liquid 
water contents present in these clouds, the 
impinging droplets do not freeze on impact and may 
run back over the accreting surface as 
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liquid water before freezing further downstream on 
the surface. This type of ice accretion is 
characterized as glaze ice and often the resulting 
ice shape displays two pronounced growth peaks, or 
horns, on either side of the stagnation line(see 
figure 1 also). The most severe aircraft 
performance degradation due to i c i ng 2 i s typically 
associated with glaze ice formations . 


Figure 2 illustrates schematically how recent 
attempts 3 " 5 to analytically model aircraft icing 
have decomposed the ice accretion process. First, 
the aerodynamic flowfield around the body of 
interest is calculated, usually by an inviscid 
panelling method. Droplet trajectories around the 
body are then calculated by integrating the 
droplet equations of motion within the body 
flowfield. From these trajectory calculations the 
mass flux of impinging droplets at a given 
location on the body can then be determined. The 
third, and crucial step, involves a thermodynamic 
analysis of the freezing process at the icing 
surface. A steady-state energy balance involving 
the impinging heat load from the droplets and the 
heat removed by convection, evaporation and 
sublimation is typically applied to small control 
volumes along the icing surface. The fourth step 
then calculates the amount of ice formed at each 
location on the body as a result of satisfying 
this energy balance and constructs the ice shape 
on the body. This entire process may then be 
repeated using the iced geometry as input for the 
flowfield calculation. Due to this form of 
time-stepped solution, any errors or inaccuracies 
in each of these four steps tend to propagate and 
may result in an unrealistic ice shape. 



Step 1 


Step 2 


Step 3 


Step 4 


Figure 2. Schematic breakdown of analytical ice 
accretion modeling process. 
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Experimental results have confirmed the 
accuracy of the flowfield and droplet trajectory 
calculations. For rime ice growth, where the 
droplets freeze at the point of impact, the 
predicted ice shapes are in good agreement with 
experimentally measured ice shapes. However, for 
glaze ice conditions the predicted ice shapes are 
extremely sensitive to the assumed convective heat 
transfer coefficient distribution over the body, 
and, to date, no analytic model has accurately 
predicted ice shapes throughout this important 
icing regime. 

Because of the sensitivity to heat transfer, 
experimental measurements of the heat transfer 
coefficient distribution over an iced surface are 
critical for validating these analytic models. 

Due to experimental difficulties associated with 
these measurements on actual ice surfaces, local 
heat transfer coefficient measurements have only 
been made around wooden or foam models of typical 
ice shapes . In addition, experimental 
measurements comparing local heat transfer 
coefficients obtained in icing wind tunnels and in 
flight are essential if natural icing conditions 
are to be accurately simulated in icing wind 
tunnels. Very little experimental data exists in 
this area. 

This paper outlines the steady-state energy 
balance for the stagnation region of a body. From 
this energy balance the condition for transition 
from rime or "dry" ice growth to glaze or "wet" 
ice growth is determined in terms of the locally 
impinging mass flux and the local heat transfer 
coefficient. By measuring the locally impinging 
mass flux and whether the accreting ice surface is 
wet or dry, it is thus possible to infer limits on 
the magnitude of the heat transfer coefficient at 
the ice surface. 

The detection of liquid water on an accreting 
ice surface is made by a unique application of an 
ultrasonic pulse-echo technique and is described. 
Results of tests employing this ultrasonic 
technique for cylinders exposed to artificial 
icing conditions in the NASA Lewis Icing Research 
Tunnel (IRT) and natural icing conditions in 
flight from the NASA Lewis Twin Otter Icing 
Research Aircraft are then presented. These 
experimental results are us?d to compare local 
heat transfer coefficients measured around a bare 
cylinder with different levels of freestream 
turbulence and surface roughness. 


STEADY-STATE THERMODYNAMIC MODEL FOR AN ICING 



The thermodynamic analysis presented in this 
paper for a surface accreting ice follows.the 
earlier work of Messinger iU and others 11 * . 

Figure 3 shows the principle modes of energy 
transfer associated with an icing surface. Heat 
is added to the surface primarily due to the 
latent heat of fusion released as the droplets 
freeze, but also from aerodynamic heating and, to 
an even smaller extent, from the kinetic energy of 
the droplets impacting the surface. Heat is 
removed from the surface primarily by convection, 
and to a lesser degree by sublimation (when the 
surface is dry) or evaporation (when the surface 
is wet). In addition, heat is absorbed from the 
surface as the supercooled droplets impinge and 
warm to 0*C. 
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Figure 3. Modes of energy transfer for an 
accreting ice surface. 


Figure 4 depicts the local control volume 
examined in this analysis. Since only the 
stagnation region of the body is considered in 
this paper, it is assumed that the only liquid 
entering the control volume is due to the 
impinging droplets. Liquid may, however, flow out 
of the control volume along the ice surface. Thus 
if the ice surface in the stagnation region is 
wet, insufficient heat is being removed to freeze 
all of the impinging liquid, and therefore the 
"freezing fraction", n, is less than unity, where 
n is given by 

M H 

n = frozen (1) 

m" 

"impinging 

When the ice surface is dry the freezing fraction 
is unity. The steady-state assumption requires 
that the rate at which energy is added to the 
control volume equals the rate at which it is 
removed, i.e. 


Q " IN " Q "0UT (2) 

At steady-state it is assumed that the ice surface 
achieves a locally uniform equilibrium 
temperature, T *. Conduction into the ice is 
assumed to be zero and chordwise conduction 
between adjacent control volume is neglected. 

With these assumptions eq. (1) may be written in 
terms of its component heat terms as 


q; 


. = m"{l, + c.(o°c - t )> 

freezing f i v $urf ;; 


aero heating “ rhV ° ° 


Q" = M"V 

droplet K.E. 


Qr 


convection 


= h(T, 


surf 


TJ 


^subl/evap = hDL s/v^ p v.surf ' Py,«) 
k 


Q IN 


1 (3) 


^OUT 


n M — un r » r 

^droplet warming w a * 



Figure 4. Stagnation region control volume used 

for steady state thermodynamic analysis. 
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Note that this equation (eq. 4) assumes all of the 
impinging mass flux, M", is frozen, thus the ice 
surface is dry and the latent heat of sublimation, 
L , is used in the mass transfer term from the 
surface. 


The ice surface will first start to become 
wet when the equilibrium temperature of the 
surface rises to 0*C. The energy balance 
appropriate for the transition between dry and wet 
growth is thus 

+ VJ + rhV^ = h{(0°C - TJ + 

2 2 ^ 

^v (p v,0»c- p v.J } + «" C w ( 5 ) 

k 

In this case the freezing fraction must be very 
si ightly -less than unity, since the surface is now 
wet, although it has been taken to be unity for 
this transitional case. The latent heat of 
vaporization, L , is however used to represent the 
evaporative cooling term. 

The local mass flux, M", of droplets 
impinging on a body is given by 

M M = fjWV w (6) 

Where 8 is the "local collection efficiency" and 
is defined as the ratio of the locally impinging 
mass flux to the freestream mass flux, i.e. 

6 = Locally Impinging Droplet Flux (7) 
p Freestream Droplet Flux 

The local collection efficiency is governed by the 
ratio of the impinging droplets inertia to the 
aerodynamic drag on them due to the freestream 
flowfield disturbance created by the body. For 
most aircraft icing regimes the local collection 
efficiency is primarily a function of the cloud 
droplet median volume diameter, MVD, and the size 


of the accreting body. Typical values for the 
local collection efficiency on the stagnation line 
of a 10cm diameter (4 M ) cylinder at a freestream 
airspeed of lOOm/s (230mph) are 0.6 for a cloud 
composed of droplets 20 microns in diameter and 
0.4 for a cloud of 12 micron droplets. Thus a 
small body moving rapidly through a cloud of large 
droplets will have a much higher local collection 
efficiency than a large body moving slowly through 
a cloud of small droplets. Using eq. 6 allows the 
critical condition for transition between wet and 
dry ice growth to be expressed as 

eWVJL f + VJ + rhV* - h{ ( 0° C - TJ + 
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It is convenient to express the convective 
heat transfer coefficient, h, in terms of the 
dimensionless Nusselt number, Nu, where 

Nu = hd /q\ 


ULTRASONIC PULSE-ECHO ANALYSIS OF AN ACCRETING ICE 

Surface 

Figure 5 illustrates the principle of the 
ultrasonic pulse-echo measurement technique 
applied for an icing surface. A small 
transducer mounted flush with the accreting 
surface emits a brief compress ional wave, the 
ultrasonic "pulse", that travels through the ice 
as shown. This wave is reflected at the ice 
surface and returns to the emitting transducer as 
an echo signal. By measuring the time elapsed 
between the emission of the pulse and the return 
of the echo, the ice thickness over the transducer 
may be calculated using the appropriate speed of 
sound for ice. By repeatedly emitting pulses the 
ice thickness can be constantly measured as ice 
accretes on the body; pulses are typically emitted 
several thousand times a second, while the time 
elapsed between pulse emission and echo return is 
on the order of microseconds for ice thicknesses 
less than 5 cm (2"). 


With d being the uniced diameter or characteristic 
dimension of the body, and k the thermal 
conductivity of air. Experimental measurements of 
local heat transfer coefficient distributions are 
often presented in terms of a power law 
relationship between the Nusselt number and the 
Reynolds number. 


Nu = ARe B = A / p V d 


Uoo 


B 


( 10 ) 


Where A and B are experimentally derived 
constants. Using eqs. 9 and 10, equation 8 may be 
rewrittenin terms of a "critical" impinging 
liquid water content, BW (eq(ll) below). 
Here the product ( BW ) .. is the wet/dry threshold 

and represents the critical local ly^impinging 
liquid water content (gramms/meter **) necessary to 
produce a wet ice surface. If BW is greater than 
this critical value the ice surface will be wet 
and the freezing fraction will be less than one, 
while if BW is less than this critical value the 
ice surface will be dry. From this equation it 
can be seen that (BW) .. depends on the assumed 
values of A and B in tne^heat transfer coefficient 
model. By measuring the locally impinging liquid 
water content, BW, and whether the resulting ice 
growth is wet or dry it is thus possible to 
compare different heat transfer models (different 
A and B constants in eq. 10) and determine which 
models best predict the observed wet or dry ice 
growth. The technique used to detect the presence 
or absence of liquid water on an accreting ice 
surface is described in the next section. 



Figure 5. Ultrasonic pulse-echo measurement 
technique. 
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In addition to obtaining the ice thickness 
from the pulse-echo transit time, the condition of 
the ice surface may also be monitored via the 
characteristics of the ultrasonic echo pattern 
received from the ice surface. Specifically it 
has been found iq that the echo patterns received 
from a dry ice surface and a wet ice surface are 
markedly different. A "single" echo, 
corresponding to the ice/air interface over the 
transducer, is received from a dry ice surface. 
However during wet ice growth the presence of 
liquid water on the ice surface creates a 
different, reflective interface, namely that 
between liquid and air. Thus during wet ice 
growth the received echo pattern contains an echo 
from the ice/water interface and further, rapidly 
varying echoes from the water/air interface. The 
presence of these echoes, which vary due to the 
distortion of the water by the impinging droplets 
and flowfield, is used to determine if the ice 
growth is wet. If only a single echo is received 
from the ice surface then the ice growth is 
determined to be dry. 


EXPERIMENTAL FACILITIES AND TEST DESCRIPTION 


Research Tunnel Tests 


The first series of tests was performed in 
the NASA Lewis Research Center Icing Research 
Tunnel. A 0.102m (4.0") diameter cylinder 
instrumented with ultrasonic transducers was 
suspended vertically from the icing tunnel roof as 
shown in figure 6. The transducers were located 
on the stagnation line of the cylinder. The echo 
signals received from these transducers were 



displayed on an oscilloscope. The oscilloscope 
screen was video-taped in order to provide a 
permanent record of the time-dependent echo 
patterns. The test procedure consisted of 
lowering the tunnel temperature to the desired 
icing cloud condition with the water spray off. 
Once the tunnel temperature had stabilized, the 
water spray system was turned on in order to 
produce a cloud of droplets of the desired size 
(median volume diameter). Typically the spray 
system was activated for a six minute period, 
after which it was turned off. Photographic and 
other measurements of the iced cylinder were then 
made and the cylinder was completely de-iced 
before the start of the next run. A total of 36 
runs for 15 different icing conditions were 
performed . 


Natural Icing Tests 


A second series of tests was performed in 
natural icing conditions using a 0.114m (4.5") 
diameter cylinder similarly instrumented with 
ultrasonic transducers on the stagnation line. An 
oscilloscope was again used to display the 
received echo signals and video-taped as before. 
The cylinder was exposed to the icing cloud via an 
experiment carrier mounted in the roof of the NASA 
Lewis Icing Research Aircraft, (De Havilland Twin 
Otter) shown in figure 8. When deployed, the 
cylinder was located 0.53n (21") into the 
freestream above the roof-line of the aircraft. 

The cylinder was exposed tnroughout the icing 
encounter, which typically lasted approximately 20 
minutes. Throughout the exposure other 
instruments mounted on the research aircraft 
measured aircraft and icing cloud parameters. In 
particular a Johnson-Wi 1 1 iams hot-wire probe was 
used to record cloud liquid water content and a 
forward scattering laser probe (FSSP) was used to 
measure the cloud droplet size distribution. Four 
separate exposures were conducted in this series 
of flight tests. 



front view 


Figure 7. Test cylinder installation for natural 

Figure 6. Test cylinder installation for icing icing flight tests on the NASA, Twin 

wind tunnel (IRT) tests. Otter, Icing Research Aircraft. 
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Experimental Results & Discussion 


The limited time available for testing in the 
icing research tunnel and the research aircraft 
prevented a precise determination of the local 
heat transfer coefficient by means of the 
experimentally measured wet or dry ice growths 
alone. However, the wet/dry ice growth data 
obtained can be used to compare the ice growth 
predicted by different heat transfer coefficients 
in the form Nu=A Re 0 and the quasi steady-state 
icing model. Recent experimental heat transfer 
coefficient measurements about a bare cylinder by 
Van Fossen et al . , in support of NASA Lewis' icing 
research programme are compared in this way. The 
bare cylinder is considered because most of the 
wet/dry measurements were obtained for relatively 
small ice thicknesses, and hence the geometry is 
essentially that of a bare cylinder. In the Van 
Fossen study heat transfer coefficients were 
measured for two different freestream turbulence 
levels, 0.5% and 3.5%, and for two different 
cylinder surface conditions, one smooth and one 
roughened with grains of sand with an average 
element height of 0.33mm. The 0.5% freestream 
turbulence level was the minimum turbulence level 
achievable in the Van Fossen tests. The 3.5% 
turbulence was chosen to characterize the higher 
turbulence level believed to exist in the icing 
research tunnel. 

Icing Research Tunnel Results 


Figure 8 shows the ultrasonically measured 
ice growth for six different icing conditions in 
the icing research tunnel. The freestream 
velocity was the same, 102.8m/s (230mph), for all 
six runs shown. The impinging liquid water 
content, BW, was determined from the 
ultrasonically measured dry ice growth accretion 
rate, since when the freezing fraction is unity 
the local ice accretion rate, d, is given by 

d ’ - S W V ~ (12) 

°ice 

In all cases the ice density, o was taken to 
be that of the pure substance. 


Figure 8 also shows the four wet/dry 
threshold curves calculated using the Van Fossen 
heat transfer coefficients. These curves are 
plotted versus ambient temperature and were 
calculated for a freestream velocity of 102.8m/s 
( 230mph) , and a cylinder diameter of 0.102m (4"). 
The four curves shown thus represent the 
transition line between wet and dry ice growth 
calculated for the four different local heat 
transfer coefficients implied by the Van Fossen 
data. If the local impinging liquid water content 
exceeds this value of ( 0 W ) .. for a given ambient 

temperature then the ice growth is calculated to 
be wet, and if BW is less than (BW) .. the ice 
growth is predicted to be dry. Fromtne figure, 
it can be seen that the heat transfer coefficient 
that best predicts the experimentally observed 
pattern of wet and dry ice growth is that measured 
for the cylinder roughened with sand and at a 
freestream turbulence level of 3.5%. While dry 
ice growth was observed at -23 * C ( - 1 0 * F ) and an. 
impinging liquid water content equal to 0.47g/m , 
the heat transfer coefficients for the 0.5% 
freestream turbulence level clearly imply wet 
growth for these conditions. Thus it appears that 
the 0.5% turbulence level heat transfer 


Cylinder Dia. = 0T02m 
Freestream Vel. = 102 8 m/s 
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Open symbols = Wet growth 



Figure 8. Plot of impinging liquid water content 
versus cloud temperature showing 
ultrasonically measured wet/dry ice 
growth and theoretical wet/dry threshold 
curves for four different heat transfer 
coefficients. (V. = 102.8m/sec 
(230mph). 

coefficients underpredict the actual heat transfer 
in the icing tunnel are therefore too low. 

Since only the stagnation region of the 
cylinder is considered, the local heat transfer 
coefficient in this region is less sensitive to 
the surface roughness than to the freestream 
turbulence level, as can be seen from the four 
curves in figure 8. However local surface 
roughness does play a critical role in determining 
where boundary layer transition occurs, and this 
in turn significantly affects the heat transfer 
distribution around the body and therefore the 
resulting ice shape. Accurate analytic models of 
the local surface roughness on real iced surfaces 
have not yet been developed; however the local 
surface roughness has been found to vary with both 
position on the body and the icing conditions 
under which the ice was *ormed. For example, a 
surface roughness of 2mm has been measured at -8*C 
(♦18*F) for ice accreted on a 0.53m chord NACA 
0012 airfoil, while at -26*C (-15*F) the surface 
roughness was found to be less than 0.1mm. Thus 
while the heat transfer coefficient implied by the 
"rough" surface, 3.5% turbulence level 
measurements appears to best predict the 
experimentally observed pattern of wet and dry ice 
growth, this heat transfer "model" may well only 
be applicable to the stagnation region and a 
different model may apply elsewhere on the 
cylinder. 
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A clear understanding of the stagnation 
region heat transfer is essential for ice shape 
prediction by current analytic icing models. 

Since these models assume that all the impinging 
liquid not frozen in the stagnation region flows 
downstream into the adjacent control volumes (see 
fig. 4), errors in the assumed stagnation region 
heat transfer coefficient will significantly 
change the mass flux distribution around the body 
during wet growth. If for example the stagnation 
region heat transfer coefficient used in such an 
icing model is less than the actual value, then 
during wet growth the model will overpredict the 
amount of liquid running out of the stagnation 
region into the adjacent control volume. This in 
turn may erroneously produce further liquid 
runback from this control volume, when in fact dry 
growth is the correct condition for this region. 
The ice shapes predicted in such a case will be in 
poor agreement with the actually observed shape, 
particularly when the feedback process inherent in 
the t i me -stepped flow field and droplet trajectory 
calculations is taken into account (fig. 2). 

Since the rough surface, 3.5% freestream 
turbulence heat transfer coefficient appeared to 
best approximate the actual heat transfer 
. occurring in the icing research tunnel (based on 
the wet/dry ice surface date from the ultrasonic 
tests), this heat transfer power law model 

Nu = 0.2460 Re 0,6444 ( 13 ) 

was compared with other ultrasonic wet/dry ice 
growth data obtained at different tunnel icing 
cloud conditions. Figure 9 shows two plots 
similar to figure 8, with impinging liquid water 
content plotted versus ambient temperature. 
Ultrasonic wet/dry ice growth measurements are 
shown for two freestream velocities, 71 . 5m/s 
(16 Omph ) and 49.2m/s (llOmph). The impinging 
liquid water content, BW, was obtained from the 
dry growth accretion rate, where available, or by 
calculation of 8 for those cases where only wet or 
marginally wet ice growth was observed. The 
wet/dry threshold curves calculated using the 3.5% 
turbulence, rough surface heat transfer 
coefficient are plotted. 

The 3.5% turbulence level, rough surface heat 
transfer coefficient correctly predicts all three 
observed dry ice growth cases. In addition, the 
two marginally wet cases observed suggest that the 
heat transfer over the ice surface may in fact be 
even greater than that implied by the 3.5% 
turbulence model. This would be consistent with 
the experimental results shown in figure 8 where 
dry growth was observed in one case while the 3.5% 
turbulence heat transfer coefficient predicted 
si ightly wet growth. 

Based on the results presented in figures 8 
and 9 it appears that the heat transfer 
coefficient model that best approximates the 
actual heat transfer occuring in the icing 
research tunnel is the 3.5% turbulence level, 
rough surface model. The actual heat transfer 
coefficient is clearly greater than those 
applicable at the low (0.5%) turbulence level and 
may be even greater the high (3.5%) turbulence 
level value. 
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Figure 9. Plot of impinging liquid water content 
versus cloud temperature showing 
ul trasonical ly measured wet/dry ice 
growth at two additional free stream 
velocities and the theoretical wet/dry 
threshold curve for the 3.5% turbulence 
level, rough surface, heat transfer 
coefficient. 
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While constant icing conditions were 
maintained throughout each exposure in the icing 
research tunnel, the natural icing cloud 
conditions, most noticeably the liquid water 
content, were not constant throughout each flight. 
Figure 10 shows a plot of cloud liquid water 
content versus exposure time for research flight 
85-24. The liquid water content was measured by a 
Johnson-Wi 11 iams hot-wire probe located near the 
nose of the aircraft (see figure 7). Also shown, 
are the experimentally observed periods of dry, 
wet and transitional ice growth, produced by the 
varying impinging liquid water content. The 
ultrasonic echo patterns received from the 
accreting ice surface were used to determine if 
the ice growth was wet, dry or transitional. Ice 
growth was characterized as transitional when the 
ultrasonic echoes received from the ice surface 
displayed neither predominantly dry nor wet echo 
patterns. 


Since the thermodynamic model used to compare 
heat transfer coefficients assumes a quasi 
steady-state icing process, and since the cloud 
liquid water content varied significantly 
throughout each exposure, a simple, steady-state 
"time constant" was estimated. This time constant 
was then applied as a criteria for comparing the 
ultrasonic data with the wet/dry growth regimes 
predicted by the steady-state model. A time 
constant of ten seconds was selected based on the 
transient thermal response of a thin ice layer and 
on the time response of the Johnson-Wi 1 1 iams 
liquid water content measurements. Therefore only 
liquid water content levels sustained for greater 
than ten seconds were used in the categor ization 
of ice growth as wet, dry or transitional. 


Figure 11 is a plot of impinging liquid water 
content versus cloud temperature. The 
experimentally observed ice growth regimes during 
the four research flights conducted are shown. 

Note that during flights 85-24 and 85-25 the full 
range of ice growth regimes was encountered with 
periods of dry, transitional and wet ice growth 
observed. No dry ice growth was observed during 
flights 85-22 and 85-23. Also shown are the four 
wet/dry threshold curves calculated using the four 
different Van Foss en heat transfer coefficients 
(0.5% turbulence, rough and smooth surface; 3.5% 
turbulence, rough and smooth surface). These 
curves were calculated for the test cylinder 
diameter of 0.114m (4.5%") and for the average 
(85-22 to 85-25) flight airspeed of 71.4m/s 
(160mph) and the average exposure altitude of 
1613m (5292'). 

Figures 10 and 11 illustrate the considerable 
variations encountered in natural icing 
conditions, both during a particular flight and 
between flights conducted on different days. For 
example, the cloud temperature, liquid water 
content and droplet size were roughly comparable 
for flights 85-24 and 85-25. However, different 
ranges of wet and dry ice growth were observed, as 
indicated by the overlapping experimental wet and 
dry growth ranges at the same impinging liquid 
water content. The implication is that the heat 
transfer differed between the two flights, both at 
nominally similar icing conditions, but conducted 
on different days through different clouds. 


Flight 85-24 
Dry = Dry ice growth 
Wet = Wet ice growth 
Tr. = Transitional ice growth 



EXPOSURE TIME, t (MINUTES) 

Figure 10. Plot of liquid water content (measured 
by the Johnson Williams probe) versus 
exposure time for flight 85-24 showing 
typical fluctuations observed in 
natural icing conditions. Also shown 
are ultrasonic measured periods of wet, 
dry, and transitional ice growth. 



Figure 11. Plot of impinging liquid water content 
versus cloud temperature showing wet, 
dry and transitional ice growth regimes 
observed in flights and theoretical 
wet/dry threshold curves for four 
different heat transfer coefficients. 
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From the figure it can be seen that for 
flight 85-25 the high (3.5%) turbulence level heat 
transfer coefficients overpredict the observed 
heat transfer, based on the steady-state model 
analysis. For this flight the low (0.5%) 
turbulence level heat transfer coefficients appear 
appropriate since both of these coefficients 
correctly predict the observed wet and dry 
regimes. The actual turbulence level applicable 
could be even less than 0.5% based on the location 
of the observed wet and dry growth regimes. 

The experimental ly observed ice growth 
regimes during flight 85-24 are consistent with 
the wet/dry threshold predicted by the 3.5% 
turbulence level, rough surface heat transfer 
coefficient. For this flight the low (0.5%) 
turbulence level models incorrectly predict wet 
growth for impinging liquid water content levels 
where dry growth was experimentally observed. 

Thus, in contrast to flight 85-25, the low 
turbulence level appears to be too low, and the 
3.5% turbulence level model gives acceptable 
results. 

The results of these tests highlight the 
effect of variations in icing conditions which are 
inherent in all natural icing encounters. Ice 
growth may thus vary from wet to dry during a 
particular encounter as the impinging liquid water 
content fluctuates. The ice shapes formed during 
pure dry (rime ice) or pure wet (glaze ice) growth 
are markedly different (see fig. 1). Therefore 
the use of an average cloud liquid water content 
to represent a particular natural icing encounter 
will in general be inappropriate since this will 
imply either pure wet or pure dry growth, when in 
reality the ice growth (wet or dry) varies as a 
function of time. 

The results of the natural icing tests also 
indicate the heat transfer occurring in natural 
icing conditions may vary from day to day despite 
similar icing conditions. One reason for this 
variation may be due to different icing cloud 
turbulence levels. Based on the limited amount of 
flight test data available it appears that in 
general the appropriate turbulence level for 
natural icing conditions is somewhat lower than 
the 3.5%+ level inferred for the icing research 
tunnel. This result is consistent with previous 
experimental comparisons between the icing 
research tunnel and flight. For these reasons, 
care should be taken in extrapolating the results 
of icing wind tunnel tests to "similar" natural 
icing cloud conditions. 

CONCLUSIONS 

Results of tests conducted in artificial 
(icing research tunnel) and natural (flight) icing 
conditions have shown the following; 

1. The presence of liquid water on an accreting 
surface may be detected using an ultrasonic 
pulse-echo technique. 

2. The threshold between dry and wet ice growth 
can be theoretically determined if the local 
convective heat transfer coefficient is known. 


3.. By comparing experimentally measured wet or dry 
ice growth with wet/dry thresholds predicted by 
different heat transfer coefficients it is 
possible to infer heat transfer coefficient 
values for an iced surface in icing conditions. 

4. The heat transfer occurring during initial ice 
growth in the icing research tunnel appears to 
be best modelled by Van Fossen's 3.5% 
freestream turbulence level, rought surface 
heat transfer model for a bare cylinder. The 
actual heat transfer coefficient may be 
slightly in excess of the values predicted by 
this model . 

5. During natural icing cloud encounters 
conditions are not constant and as a result 
periods of wet, dry and transitional ice growth 
may be observed within a single encounter. 

6. The heat transfer occuring during initial ice 
growth in natural icing conditions has been 
inferred to vary between that predicted at the 
3.5% Freestream turbulence level and the 0.5% 
level, using Van Fossen's data. 

7. Due to variations in natural icing cloud 
conditions care should be taken in 
extrapolating results from icing wind tunnels 
to "similar" natural icing cloud conditions. 
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